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We present a simple novel construction of recursion operators for integrable multidimensional 
dispersionless systems that admit Lax-type representations whose operators are linear in the spec¬ 
tral parameter and do not involve the derivatives with respect to the latter. New examples of 
recursion operators obtained using our technique include inter alia those for the general heavenly 
equation, which describes a class of anti-self-dual solutions of the vacuum Einstein equations, and 
a six-dimensional equation resulting from a system of Ferapontov and Khusnutdinova. 


Introduction 

The field of integrable systems is a dynamic subject closely related to many parts of modern mathematics 
and featuring manifold applications, cf. e.g. mK00@liai2DlBEll3Z10DlHZ] and references therein. 

A characteristic feature of integrable systems is that they never show up alone - instead, they always 
belong to integrable hierarchies, and the latter are typically constructed using the recursion operators 
(ROs), cf. e.g. [T1 IHl [T2[ I2T113611391140] and references therein. In particular, the ROs for multidimensional 
integrable dispersionless systems are now a subject of intense research, cf. e.g. m H2J EU EHl [3011331 EH [351 
[40] and references therein. In fact, the overwhelming majority of integrable systems in four or more inde¬ 
pendent variable known to date, including those relevant for applications, like the (anti-)self-dual Yang- 
Mills equations or the (anti-)self-dual vacuum Einstein equations, is dispersionless, i.e., can be rewritten 
as first-order homegeneous quasilinear systems, cf. e.g. PUZl [45] and Example 5 below for details. 

The study of a hierarchy associated to a given integrable system, and hence of its recursion operator, 
is important for many reasons, inter alia because known exact solutions of integrable systems, like, say, 
multisoliton and finite-gap solutions, always happen to be invariant under a certain combination of the 
flows in the hierarchy, and this is instrumental in the actual construction of the solutions in questions, 
cf. e.g. [121136]. 

Another feature making recursion operators relevant for applications is that existence of a (formal) 
recursion operator can in many situations be employed as an efficient integrability test or as a means 
of proving nonintegrability, see e.g. [32] and references therein for general theory and [4S.J for a recent 
application. 

There are several methods for construction of the ROs for the integrable dispersionless systems: the 
partner symmetry approach [24], the method based on adjoint action of the Lax operators [3U] 5 and the 
approach using the Cartan equivalence method [331 [341 135] . 

Below we present another method for construction of ROs for integrable multidimensional disper- 
sionlcss systems. Roughly speaking, it consists in constructing a Lax-type representation for the system 
under study with the following additional properties: i) the Lax operators are linear in the spectral 
parameter and ii) the solution of the associated linear system is a (nonlocal) symmetry of the nonlinear 
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system under study. Taking the coefficients at the powers of spectral parameter in the operators con¬ 
stituting the newly constructed Lax-type representation yields, under certain technical conditions, the 
RO. A similar approach, with symmetries replaced by cosymmetries and Proposition [Tj by Proposition [2l 
allows one to look for the adjoint ROs. 

To the best of author’s knowledge, the key idea of our approach, i.e., that the Lax-type representation 
employed for the construction of the recursion operator should be constructed from the original Lax-type 
representation but is by no means obliged to be identical with the latter, has not yet appeared in the 
literature. This kind of flexibility noticeably enhances the applicability of our method. On the other 
hand, the extraction of the RO from a Lax-type representation which is linear in the spectral parameter 
was already explored to some extent e.g. in ra. 

The paper is organized as follows. In Section [Tj we recall some basic facts concerning the geometric 
theory of PDEs. In Section [2] we state some general results on the construction of recursion operators. 
Section 3, which is the core of the paper, explores the construction of special Lax-type representations 
that give rise to ROs through the results of Section [21 Section 0 ] gives a selection of examples illustrating 
our approach, and Section O provides a brief discussion. 


1 Preliminaries 


Let F = 0 be a system of m PDEs Fj — 0, / = 1,..., m, in d independent variables x l , i — 1,..., d, so 
we put x = (x 1 ,... ,x d ), for an unknown iV-component vector function u = (it 1 ,... ,u N ) T , where the 
superscript ‘T’ indicates the transposed matrix. Let 

ul in = a il+ -+^n“/a(x 1 ) il • • • d(x n ) in 


and Moo - o = I 11 what follows all functions are tacitly assumed to be meromorphic in all their 

arguments. 

As e.g. in uni H3, x l and uf A are considered here and below as independent quantities and can 
be viewed as coordinates on an abstract infinite-dimensional space (a jet space). By a local function we 
shall mean a function of x, u and of a finite number of the derivatives of u. To indicate locality of a 
function / we shall write it as f(x, [u]). 

We denote by 


N 


D X ] — 


dxi 


d i v oo <-\ 

Af V v° d 

/ j / j u, i 1 ...ij- 1 (ij+l)ij +1 ...i n Q u a 
a=l *i,...,* n =0 h-.-in 


the total derivatives with respect to ad, cf. e.g. mm For the sake of simplicity the extensions of D x <. 
to nonlocal variables are again denoted by D x% , as e.g. in [ 3U] . 

The system F = 0 together with all its differential consequences 


Df, ■ ■ ■ D l " n T = 0, zi,i 2 , • • • An = 1, • • • ,d, n — 1,2, , 

defines a diffiety Soljr, a submanifold of the infinite jet space, see e.g. m and references therein for 
details. This submanifold can be informally thought of as a set of all formal solutions of F = 0, which 
motivates the choice of notation. Below all equations will be required to hold only on Soljr, or a certain 
differential covering thereof, see e.g. 031 for details on the latter, rather than e.g. on the whole jet space 
unless otherwise explicitly stated. 

The directional derivative (cf. e.g. [3, [30] and references therein) along U = (U 1 ,..., U N ) T is the 
vector field on the jet space 


N OO 

sv = Y. E F 

a= 1 


■D^ n U a ) 


d 


dU t ...in 
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The total derivatives as well as the directional derivative can be applied to (possibly vector or matrix) 
local functions P. 

Recall, cf. e.g. mm , that a local TV-component vector function U is a (characteristic of a) symmetry 
for the system T = 0 if U satisfies the linearized version of this system, namely, £ F {U) = 0. We stress 
that by our blanket assumption this is required to hold on SoIjf only rather than on the whole jet space. 
Informally this definition means that the flow u T = U leaves Solj- invariant. 

Here 



N OO 


EE 


8{ 


dU°: 


1 1 • • •'In 


D\ ■■■D 

X L 


in 

X n 


is the operator of linearization and £j = (i Fl ,... ,£ Frn ) T . Notice an important identity £f(U) = djj(f) 
which holds for any local / and U, see e.g. Ha- 

Finally, recall that if we have an operator in total derivatives 


Q 


q d d 

q°+ EE'-'E^ 1 '*- 0 ^"'- 0 ^ 

fc=i ji=i j k =i 


where Q° and Q n '" 3k are s x .V-matrix-valued local functions, its formal adjoint reads (cf. e.g. [4U]) 


Q f = (q«) t +e E - ■ Ef- 1 ) 1 - 0 - ■ ■ ■ a* ° 

k =1 j i=l jk =1 


For integrable systems in more than two independent variables their symmetries usually depend, in 
addition to local variables, described above, also on nonlocal variables. In other words, we are forced 
to consider solutions U of £ F (U) = 0 involving nonlocal variables. A number of authors refers to such 
objects as to the shadows of nonlocal symmetries, see e.g. mi and references therein for more details 
on this terminology; cf. also [5] and references therein. However, as in the present paper we shall not 
deal with full nonlocal symmetries in the sense of mi. for the sake of simplicity in what follows we shall 
refer to solutions U of t F {U) = 0 involving nonlocal variables just as to nonlocal symmetries. 

A cosymmetry , also known as adjoint symmetry, see e.g. [Sj, 7 is a quantity which is dual to a 
symmetry: it has m rather than N components and satisfies (cf. e.g. [31 |T?j for details) the system 

h( 7) = 0. 


Note that cosymmetries, just like symmetries, may depend on nonlocal variables in addition to the local 
ones. 

A recursion operator (RO) for our system T = 0 is then, roughly speaking, an operator that sends 
any symmetry of T = 0 into another symmetry thereof [3HJ 0Uj. However, this definition works well 
only when the RO contains no nonlocal terms, as often is the case for linear and lincarizable equations, 
cf. e.g. mm and references therein. 

In order to properly handle nonlocal terms in ROs it is more appropriate to view the RO as a 
Backlund auto-transformation for £j^(G) = 0, see [30) E, 32j and references therein for details, and 
below we adhere to this point of view. Likewise, an adjoint recursion operator is from this perspective 
a Backlund auto-transformation for ^( 7 ) = 0, cf. e.g. [T7] . 

As a closing remark, recall that a partial differential system is called dispersionless (cf. e.g. [28j and 
references therein) if it can be written in the form of a quasilinear homogeneous first-order system 


d n 

EE-V 

j =1 « =1 


du a 


= 0 , 


N\T 


where I — 1,..., m, m > N, and u = (u,..., u N ) 


( 1 ) 
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This class of systems is quite rich: for instance, quasilinear scalar second-order PDEs which do not 
explicitly involve the dependent variable u, 

d d 

EE P 1 (x, du/dx 1 ,..., du/dx d )d 2 u/dx l dx j = 0, 

i= 1 j=i 

can be brought into the form dH) by putting u = (du/dx 1 ,..., du/dx d ) T . 


2 Some remarks on recursion operators 

Consider the following differential operators in total derivatives 

A A ~ ypn,i. r>, Ip- V i = 1,2, 

3 = 1 J =1 

d d 

L = L° + £ L k D x k, M = M° + E M k D x k, 

k= 1 k =1 


( 2 ) 


where Aj = A\(x, [it]) and Bj = B/(x, [it]) for i = 1,2 and j = 1 ,d are scalar functions, A ° = 
A/(x, [it]) and B) = B/(x, [it]) are N x N matrices, L° = L°(x, [it]) and L k = L k (x, [it]) for k = 1,..., d 
are N x m matrices, while M° = M°(x, [it]) and M k = M k (x, [it]) for k — 1,..., d are m x N matrices. 

Proposition 1. Let the operators Aj, Bi, L, M of the form (Q|) be such that 

i)[A 1 ,A 2 ]=0, (3) 

n) (A 1 B 2 -A 2 B 1 ) = L oi T , (4) 

Hi) ijr = M o (BiA 2 — B 2 Ai), (5) 

iv) there exist two distinct indices p, q G { 1 ,..., d} such that we can express D xP U and D x qU from the 
relations 

A l (U) = B l (U), t = l, 2 . ( 6 ) 

Then relations (OJ) define a recursion operator for T = 0, i.e., whenever U = (U 1 ,..., U N ) T is a 
(possibly nonlocal) symmetry for T = 0, so is U = (U 1 , ..., U N ) T defined by @). 

Proof. First of all, if U — (U 1 ,..., U N ) T is a symmetry for T = 0, then the system ([ 6 ]) for U is 
compatible by virtue of i)-iii). 

We now have £jr(U) = 0 by virtue of ([5]), and hence U is a shadow of symmetry for our system. □ 
As an aside note that the condition (J5D in a slightly different form has appeared in |TT, T5] [TB] , and 
was given there a certain geometric interpretation. 

In complete analogy with the above we can also readily prove the counterpart of Proposition [1] for 
adjoint recursion operators. 

Proposition 2. Suppose that Ai, B iy L , M are as before, but with A® and B/ being m x m matrices. 
Further assume that these operators are such that 

i)[A u A 2 ]= 0 , (7) 

H) P T = Lo(B x A 2 -B 2 Afi), ( 8 ) 

Hi) (AiB 2 - A 2 BP = M O 4, (9) 

iv) there exist two distinct indices p, q e { 1 ,..., d} such that we can express D xP 7 and D x >from the 
system 

A i (p) = B i ( 1 ), i — 1,2. (10) 

Then m defines an adjoint recursion operator for J 7 = 0 , i.e., whenever 7 = ( 7 1 ,... , r ) m ) T is a 
cosymmetry for T = 0 , then so is 7 defined by ITTh) . 
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3 Recursion operators from Lax-type representations 

We start with the following observation which is readily checked by straightforward computation. 

Proposition 3. Suppose that the operators Ai and Bi define a RO (resp. an adjoint RO) as in Propo¬ 
sition 0 (resp. as in Proposition dj). Then the operators £, = A Ai — Bi, i = 1,2, where X is a spectral 
parameter, satisfy [C\, Cf\ = 0, i.e., they constitute a Lax-type representation for T = 0. 

Thus, if our system admits an (adjoint) RO, it also admits a Lax-type representation which is linear 
in the spectral parameter. 

Hence, a natural source of Ai, B, and L, M satisfying the conditions of Proposition [Tj is provided by 
the Lax-type representations for T = 0 of the form 

Ciif = 0 , * = 1 , 2 , ( 11 ) 

with Ci linear in A such that if is a nonlocal symmetry of T = 0, i.e., we have MVO = °> cf - e -S- 
[4lj. Then putting Ci = XBj — Ai or C % = A Ai — B, gives us natural candidates for Ai and Bi which 
then should be checked against the conditions of Proposition [Q and, if the latter hold, yield a recursion 
operator for T = 0. 

Even if it is impossible to construct the recursion operator from (TIT]) , one still can construct infinite 
series of nonlocal symmetries for T = 0 expanding if into formal Taylor or Laurent series in A, cf. e.g. 

mm- 

It is important to stress that the Lax representation with the operators C % employed in the above con¬ 
struction does not have to be the original Lax representation of our system. In general, we should custom 
tailor the operators £ 12 constituting the Lax pair for the construction in question, so that the solutions 
of the associated linear problem (TTT1) are symmetries, i.e., satisfy the linearized version of our system. 

The natural building blocks for these C, are the original Lax operators SPi, their formal adjoints 
PP-f , and the operators ad.^r. = \SPi, •], but in general one has to twist them, cf. e.g. Example 4 below. 

More explicitly, suppose that the system under study, i.e., T = 0, admits a Lax representation of 
the form 

SPpf = 0, * = 1,2, (12) 

where SPi are linear in the spectral parameter A but such if does not satisfy £p{if) = 0, i.e., if is not a 
(nonlocal) symmetry. 

Then we can seek for a nonlocal symmetry of T = 0 of the form 

$ = $(£, [u,if,x,(])i ( 13 ) 

i.e., for a vector function of x, and of u , if, and ( and a finite number of the derivatives of u, if, y, 
and <f. 

Here y satisfies the system 

SX'iX = 0 , * = 1 , 2 , 

d 

% = V (W/dxi satisfies 

3 = 1 

['%„ %‘\ = o, * = 1,2, 

and C = (CV--,C d ) T - 

Moreover, we should also require that there exist the operators which are linear in A and such 
that 

j^$ = 0, * = 1,2. 

Then one should extract Ai and Bi for Proposition [T] from these ACi rather than from the original SPi. 
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The study of specific examples strongly suggests that the nonlocal variables ip, x and C usually do 
not mix and enter <f> linearly, so it typically suffices to look for nonlocal symmetries $ of T = 0 in one 
of the following forms, 

r r d 

<r = + Y,'E ak ‘ D ‘ t ^’ < 14 > 

s=l s =1 k= 1 

r r d 

*“ = x>'*x.+< 15 ) 

s=l s=l k =1 

d d 

(c‘), (is) 

3=1 k,l =1 

where a = 1,..., N, p G N, a“, a" ,fc , 6“ ,s , & a,s,A: , c“ and c“’ fc are local functions, and r is the number 
of components of ip, instead of general form ffl3]l . Moreover, in many cases it is possible to restrict 
oneself to considering only zero-order terms in the nonlocal variables, i.e., use the simpler Ansatze 

<f) Q = agip s , <3>“ = b a,s Xs , and c< jC j instead of (Thfll . (TT5T) and (TT6T) . 

s=l s=l j=l 

A similar approach, with symmetries replaced by cosymmetries and Proposition [Q by Proposition [2] 
can, of course, be applied to the construction of adjoint recursion operators. 

Finally, let us point out that the assumption of linear dependence of the Lax operators on A is 
not as restrictive as it seems: a great many of known today multidimensional dispersionless hierarchies 
include systems with this property and, moreover, by Proposition [3] if a system admits an RO of the 
form described in Proposition [T] then it necessarily possesses a Lax-type representation whose operators 
are linear in A. 

Moreover, in some cases the nonlinearity in A can be removed upon a proper rewriting of the Lax 
pair. 

Consider for example the Pavlov equation |43j 

'Uyy T U xt T Hx'U'xy ^y^xx 0, (IT) 

which possesses a Lax pair of the form [53] (cf. also [26]) 

Ipy = (-Az - tpt = (“A 2 - A U x + U y )lp x (18) 

which is quadratic in A. 

However, rewriting the second equation of (fT8|l as 

ipt = A(-A - u x )ip x + u y ip x 
and substituting ip y for (—A — u x )ip x , we obtain 

1p t = XlPy A UylPx, 

and thus (TT8|) can be rewritten in the form linear in A: 

ply = ( -U x - A )lpx, Ipt = -Xlpy + Uyipx- (19) 


4 Examples 

Example 1 

It is readily checked that if ip satisfies (1191) then $ = l/ip x is a nonlocal symmetry for (1X71) . i.e., U = l/ip x 
satisfies the linearized version of (ITT)) . 

bXyy T U xt T U J X.U xy ^yU xx -)- U xy U x U XX Uy 0. (20) 
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Now by virtue of (fl9|) $ = l/ip x also satisfies a pair of linear equations of the form — 0, % — 1, 2, 
where 

Dy T (A U X )D X T U X X1 ^2 T XDy UyD X T U X y. 

Now, as <f> = l/Vh satisfies (12UD we can identify ^ with from (HIT) and put Jd) = XAi — Bi, where 

A\ D X i A 2 Dyi U I "f~ U X B X Uxxi -^2 Dt T UyD x U X y. 

Then it is easily seen that these operators satisfy the conditions of Proposition [Q for suitably chosen L 
and M, so we arrive at the recursion operator for the Pavlov equation given by the formulas 

U x Uy T u x U x u xx U ", Uy Ui T UyU x u X yU , 

which maps a (possibly nonlocal) symmetry U to a new nonlocal symmetry U. This is nothing but the 
recursion operator found in [2SJ rewritten as a Backlund auto-transformation for ( 12 Ujh 


Example 2 

Consider the general heavenly equation 


UU X yU z t T bu X zUyt + CU x f,Uy Z 0 , (X b C 0 , ( 21 ) 

where a, b, c are constants. Here d = 3, N = 1, x 1 = x, x 2 = y, x 3 = z, x 4 = t, u l = u. Note that this 
equation describes inter alia a class of anti-self-dual solutions of the Einstein field equations as shown 
in [22] - 

By definition, a (nonlocal) symmetry U of (Ell) satisfies the linearized version of (I2T|) . that is, 


UU X yU z t T UU z tU X y T bUxzUyt -(- bzty^U xz T CUxt.Uyz T CUy Z U x t d. 


( 22 ) 


Eq. (12TT) admits [TOl E2| a Lax-type representation with the operators 


L\ 

L 2 


(1 + c\)D x 
(1 - b\)Dy 


—D t - c\—D z , 

Uzt U zt 


^D t + bX^D z . 

Uzt Uzt 


(23) 


It is readily seen that ^ that satisfies — 0, i — 1, 2, also satisfies (l22]h i.e., it is a nonlocal symmetry 
for (12B . so in spirit of Proposition [3] let 


Hi = cD x - c^P 2 , H 2 = -bD y + b^D z , 

Uzt U zt 

B 1 = —D x + —D t , B 2 = -D y + ^lD t . 

Uzt Uzt 

Then all conditions of Proposition Q] are readily seen to be satished, e.g. we have B 1 A 2 — B 2 A 1 = 
(1 /u z t)^F, where F = au xy u z t + bu xz u yt + cu xt u yz . 

Hence the relations 


U x = 


'LL XZ Uf 4 “ Cli x tU z u z tU x 


> Uy = 


Uy Z U t 


blLytUZ 


'LLztUy 


bu 


zt 


(24) 


where U is a (possibly nonlocal) symmetry of (l 2 Tlh define a new symmetry U and thus a novel recursion 
operator ()24l) for () 2 Tlh i.e., a Backlund auto-transformation for (1221) . This operator has hrst appeared in 
the very hrst version of the present paper (arXiv:1501.01955vl) and was later rediscovered, in a slightly 
different form, in [23J. 
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Example 3 

Consider the Martinez Alonso-Shabat [29] equation 

Uyt. 'U'z'U'xy 'Uy'U'xz (25) 

where d = 4, N = 1, x 1 = x, x 2 = y, x 3 = z, x 4 = t, u 1 = u. Eg. (121)1 admits pS] a Lax-type 
representation with the operators 

L± = D y — \u y D x , L 2 = D z — Xu z D x + XDf. (26) 

Consider the action of operators 

d-'i \Lii •] 

on operators of the form wD x and spell it out as 

ad Li (wD x ) = ((A Bi - Ai)w)D x , 


so we have 

A\ Dy, A 2 Dzi A1 ztyD x u x y : B2 u z D x ID ^ u 

All conditions of Proposition [Tj are again satisfied, so the relations 

Uy Uy U x U X yU , 

U z U Z U X U X zU i A ■ 


(27) 


where U is any (possibly nonlocal) symmetry of (l25]h dehne a new symmetry U, i.e., Eq. (l27]l provides 
a recursion operator for (1251) . This operator was already found in [35]. 

Applying (l27j) to the simplest symmetry U = u x we obtain (modulo an arbitrary function of x and 
t resulting from the integration) a nonlocal symmetry of the form U = w — u 2 x / 2 , where the nonlocal 
variable w is defined by the relations 


Wy r U j y'U j xxi 'll) z V'z^xx ^ xt • 

This example suggests that Eq. (l27l) can be further simplified: a new symmetry U can be constructed 
by putting 

U = W - u x U , 

where W is defined by the relations 

fPy UyU X T U X Uy, 

W z u z U x T u x U z U t . 

Thus, the recursion operators produced within our approach are not necessarily in the simplest possible 
form. 


Example 4 

Consider a system da nu 


m t = n x + nm r — mn r , n z = rn y + mn s — nm s , (28) 

where x, y, z, r, s, t are independent and m, n are dependent variables. 

Eq. (j28l) can be written as a condition of commutativity for the following pair of vector fields, 


[D z — mD s — XD x + XmD r , D y — nD s — XD t + XnD r \ = 0 . 


( 29 ) 


where A is a spectral parameter, and is therefore integrable. 

By virtue of ([28]) there exists a potential u such that 

m = u z /u s , n = Uy/u s . (30) 

Substituting this into (1281) gives a single second-order equation for u, 


U s U z t U z U s t U s U X y “t“ 'Uy'U'sx l U"y'U / rz T* U z U r y 0, 


(31) 


which can be written as a commutativity condition [Li,I/ 2 ] = 0 for 

L 1 =D Z - —D s - A D x + X—D r , L-2 = D y — — D s - XD t + X—D r . 

u s u s u s u s 

Let x satisfy L\x — 0, i — 1, 2. Then Q = u s /x is a nonlocal symmetry for (T3TT) and we can readily 
obtain a linear system for £ of the form C\C, = 0, £ 2 C = 0 from that for y. 

Upon spelling out L t as L{ = A z — XBi and checking that the conditions of Theorem [[] are satisfied 
for suitable L and M we arrive at the following novel RO for (j3T[) : 


u y = ^u s - ^ U r + U t - (Mst Ury) U , 


lie 


Us 


Uc 


Uy - Uy (u r y — U o~) 

U z = —U s - —U r + U X + ALL - Liu. 

Ue Uo Uo 


(32) 


ffere again U is a (in general nonlocal) symmetry of (1311) . and U is a new nonlocal symmetry for (j3li) . 


Example 5 

Consider (see e.g. BE SB and references therein) the following Lax operators: 

L x = D y + d x K — A D x , L 2 = D x + dyK — X Dy. (33) 

The commutativity condition [L 1 ; L 2 ] = 0 yields the system 

d x d x K - dydyK - [d x K, dyK } = 0. (34) 

ffere K takes values in a (matrix) Lie algebra g and is known as the Yang Ji-matrix. 

Note that (1341) can be rewritten in the dispersionless form of the type ([I]) upon introducing (cf. e.g. 
mm) an additional dependent variable J taking values in G = exp(g). Namely, the dispersionless 
form of (j34l) reads 


d x J — JdyK , d y J = Jd x K. 

The compatibility condition d x d y J = d y d x J yields an equation 

(35) 

d x (J~ 1 d x J) = dy(J~ l d y J), 

(36) 


while (1341) is identically satisfied by virtue of (1351) . In other words, (l35j) defines a Backlund transforma¬ 
tion relating (134|) and (136|) . 

If we define the gauge potentials by setting A x = dyK , A y = d x K , A x = 0, Ay = 0, then they 
satisfy the anti-self-dual Yang-Mills equations on K 4 with Euclidean metric with coordinates X 1 such 
that \/2x = X 1 + LY 2 , \[2x = A " 1 — LA 2 , y/2y = X 3 — LA 4 , y/2y = X 3 + LA 4 , where i = >/—1, cf. e.g. 
[HI]. Conversely, any sufficiently smooth solution of the anti-self-dual Yang-Mills equations can, up to 
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a suitable gauge transformation, be obtained in this way [3Tj . In other words, the anti-self-dnal Yang- 
Mills equations, which play an important role in modern physics, cf. e.g. [31] and references therein, 
can, in a suitable gauge, be written in dispersionless form, exactly as claimed in Introduction. 

It is natural to assume that in the Lax pair equations Lj-0 = 0 we have -0 e G = exp(g). However, 
such -0 cannot be a symmetry for (154]h since a symmetry must live in g just like K. However, we can 
get a Lax representation for e g using the adjoint action: [L i, <f>] = [L 2 , $] = 0, i.e., we use the Lax 
operators Ci = ad Lj. 

Then all conditions of Proposition [T] are satisfied for Ai and Bi obtained by spelling out Ci as 
Ci = A Ai — Bi, which yields the following RO for (l34lh 

V x = U y + [d x K, U], Vy = U x + [t %K , U\ , (37) 

where U and V are symmetries for (l34lh This RO is known from the literature, see [6J 02) [2]. 

If in the above construction we take g = hiff(M iV ), i.e., we put 

N 

K = y^'Ujd/dz 1 , 

i =1 

where Ui = Ui(x, y, x, y, z 1 ,..., z N ) are scalar functions, then (134)1 becomes [28] the Manakov-Santini 
system [25], an integrable system in (N + 4) independent variables, and (1371) provides a recursion 
operator for this system which was found earlier in [3D] (see also [28j) by a different method. 

Further examples of recursion operators obtained using our approach can be found in [ID]. 


5 Concluding remarks 

We presented above a method of construction for recursion operators and adjoint recursion operators for 
a broad class of multidimensional integrable systems that can be written as commutativity conditions 
for a pair of vector fields, or, even more broadly, of linear combinations of vector fields with zero- 
order matrix operators, linear in the spectral parameter and free of derivatives in the latter. This 
method is based on constructing a special Lax-type representation for the system under study from the 
original Lax-type representation, see Section [3] for details. Note that if a system under study admits 
several essentially distinct Lax representations then our approach applied to these Lax representations 
in principle can give rise to several essentially distinct recursion operators. 

To the best of author’s knowledge, the method in question works for all known examples of mul¬ 
tidimensional dispersionless integrable systems admitting recursion operators in the form of Backhand 
auto-transformation for linearized systems, including e.g. the ABC equation [45] 

dU X Uyi -(- bUytlxi T CUiU X y = 0, fl Y i Y c = 0, 

the simplest (2+l)-dimensional equation of the so-called universal hierarchy of Martinez Alonso and 
Shabat [29], 

^yy ^y^tx Y W T 11./y 0, 

the complex Monge-Ampere equation (see [35] for its recursion operator), first, second and modified 
heavenly equations, etc. In particular, this method enabled us to find hitherto unknown recursion 
operators for the general heavenly equation of Doubrov and Ferapontov (I2TT) and for the equation 
(Bill. Note that our approach, when applicable, is in general computationally less demanding than the 
methods of [30] and [341] . 

While all examples to which we have applied our approach so far are dispersionless, Propositions [T| 
and [2] and the method described in Section [3] do not explicitly require this to be the case, so it would 
be interesting to find out whether it is possible to extend our approach to the construction of recursion 
operators to other classes of integrable systems. 
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